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Abstract 

Given a graph G, let G a be an oriented graph of G with the orientation a and skew- 
adjacency matrix S(G a ). The skew energy of the oriented graph G" 7 , denoted by £g(G a ), 
is defined as the sum of the absolute values of all the eigenvalues of S(G cr ). In this paper, 
we study the skew energy of random oriented graphs and formulate an exact estimate of the 
skew energy for almost all oriented graphs by generalizing Wigner's semicircle law. More- 
over, we consider the skew energy of random regular oriented graphs G a n d , and get an exact 
estimate of the skew energy for almost all regular oriented graphs. 
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1 Introduction 

Let G be a simple undirected graph with vertex set V(G) = {v±,V2, ■ ■ ■ ,v n }, and let G a be 
an oriented graph of G with the orientation a, which assigns to each edge of G a direction so 
that the induced graph G a becomes a directed graph. The skew- adjacency matrix of G a is the 
n x n matrix S(G a ) = [sij], where = 1 and = — 1 if (vi,Vj) is an arc of G a , otherwise 
Sij = Sji = 0. The skew energy [l] of G a is defined as the sum of the absolute values of all the 
eigenvalues of S(G a ), denoted by £s(G a ). Obviously, S(G a ) is a skew-symmetric matrix, and 
thus all the eigenvalues are purely imaginary numbers. 

Since the concept of the energy of simple undirected graphs was introduced by Gutman in [S] , 
there have been lots of research papers on this topic. We refer the survey [9] and the book [12] 
to the reader for details. The energy of a graph has a close link to chemistry. An important 
quantum-chemical characteristic of a conjugated molecule is its total 7r-electron energy. There 



'Supported by NSFC and the "973" project. 



1 



are situations when chemists use digraphs rather than graphs. One such situation is when 
vertices represent distinct chemical species and arcs represent the direction in which a particular 
reaction takes place between the two corresponding species. It is reasonable to expect that the 
skew energy has similar applications as energy in chemistry. 

Adiga et. al. [lj first defined the skew energy of an oriented graph, and obtained some 
properties of the skew energy. They derived an upper bound for the skew energy and constructed 
a family of infinitely many oriented graphs attaining the maximum. They also proved that the 
skew energy of an oriented tree is independent of its orientation, and equal to the energy of 
its underlying tree. Then, Shader et. al. [15] studied the relationship between the spectra of 
a graph G and the skew-spectra of an oriented graph G u of G, which would be helpful to the 
study of the relationship between the energy of G and the skew energy of G a . Hou and Lei [TU] 
characterized the coefficients of the characteristic polynomial of the skew-adjacency matrix of 
an oriented graph. Moreover, other bounds and the extremal graphs of some classes of oriented 
graphs have been established. In [11] and [16], Hou et. al. determined the oriented unicyclic 
graphs with minimal and maximal skew energy and the oriented bicyclic graphs with minimal 
and maximal shew energy, respectively. Gong and Xu [7] characterized the 3-regular oriented 
graphs with optimum skew energy. 

It is well known that it is rather hard to compute the eigenvalues for a large matrix and by the 
extremal graphs we can hardly see the major behavior of the invariant £s(G) for most oriented 
graphs with respect to other graph parameters. Therefore, in this paper, we will study the skew 
energy in the setting of random oriented graphs. We first formulate an exact estimate of the 
skew energy for almost all oriented graphs by generalizing Wigner's semicircle law. Moreover, 
we investigate the skew energy of random regular oriented graphs, and also obtain an exact 
estimate. 

Various energies of random graphs have been studied. Du et. al. considered the Laplacian 
energy in [3] and the energy in [6J. Moreover, they also investigated other energies in [5], such as 
the signless Laplacian energy, incidence energy, distance energy and the Laplacian-energy like 
invariant. It is worth to point out that their results depend on the limiting spectral distribution 
of a random real symmetric matrix. But our results on the skew energy of a random oriented 
graph relies on the limiting spectral distribution of a random complex Hermitian matrix. 

The rest of the paper is organized as follows: In Section 2, we will list some notations and 
collect a few auxiliary results. Then in Section 3, we will consider the skew energy of random 
oriented graphs. Finally, in Section 4, we will be devoted to estimating the skew energy of 
random regular oriented graphs. 
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2 Preliminaries 



In this section, we state some notations and collect a few results that will be used in the 
sequel of the paper. 

Given an Hermitian matrix M on order n, denote its n eigenvalues by 

X 1 (M),X 2 (M),...,X n (M), 
and the empirical spectral distribution (ESD) of the matrix M by 

Fn{x) = - {\i(M)\\(M) <x,i = l,2,...,n} 

n 

where 1 1\ means the cardinality of the set /. The distribution to which the ESD of the random 
matrix M converges as n — > oo is called the limiting spectral distribution (LSD) of M. 

The study on the spectral distribution of random matrices plays a critical role in estimating 
the skew energy of random oriented graphs. One pioneer work in the field of the spectral 
distribution of random matrices [2 jll4j is Wigner's semicircle law discovered by Wigner in |18pi9| . 
which characterizes the limiting spectral distribution of a sort of random matrices. This sort of 
random matrices is so-called the Wigner matrices, denoted by X n , which satisfies that 

• X n is a Hermitian matrix, 

1 < i < j < n , where Xji means the conjugate of 

Xjij 

• the upper-triangular entries Xij, 1 < % < j < n, are independently identically distributed 
(i.i.d.) complex random variables with mean zero and unit variance, 

• the diagonal entries xu, 1 < i < n, are i.i.d. real random variances, independent of the 
upper-triangular entries, with mean zero, 



oo. 



• for each positive integer k, max (E(|xn| fc ), E(|a^i2 1 fe )) < 

Then the Wigner's semicircle law can be stated as follows. 

Theorem 2.1 ]19$ Let X n be a Wigner matrix. Then the empirical spectral distribution 
F™ 1/2Xn (x) converges to the standard semicircle distribution whose density is given by 



p sc (x) := -^-y/A- x 2 



2vr v - - M<2- 

Given a random graph model G(n,p), we say that almost every graph G(n,p) £ Q{n,p) has 
a certain property V if the probability that G(n,p) has the property V tends to 1 as n — > oo, 
or we say G(n,p) almost surely (a.s.) satisfies the property V. In the sequel, we shall consider 
two random graph models: the random oriented graph model Q a (n,p) and the random regular 
oriented graph model G^ d , the definitions of which will be given later. 
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In this paper, we use the following standard asymptotic notations: as n — > oo, f(n) = o(g(n)) 
means that f(n)/g(n) — > 0; /(n) = ui{g{n)) means that f(n)/g(n) — > oo; /(n) = 0{g(n)) means 
that there exists a constant C such that |/(n)| < Cg(n); f(n) = Q(g(n)) means that there exists 
a constant c > such that /(n) > cg(n). 



3 The skew energy of G a {n,p) 

In this section, we consider random oriented graphs and obtain an estimate of the skew 
energy for almost all oriented graphs by generalizing Wigner's semicircle law. 

We first give the definition of a random oriented graph G a (n,p). Given p = p(n), < p < 1, 
a random oriented graph on n vertices is obtained by drawing an edge between each pair of 
vertices, randomly and independently, with probability p and then orienting each existing edge, 
randomly and independently, with probability 1/2. That is to say, for a given oriented graph 
G = G a (n,p) with m arcs, P(G) = p m (l — p)(^~ m ■ 2~ m . Apparently, the skew- adjacency 
matrix S(G a (n,p)) = [sij] (or S n , for brevity) of G a (n,p) is a random matrix such that 

• S n is skew-symmetric, i.e., = —Sji for 1 < i < j < n, and in particular, su = for 

1 < i < n; 

• the upper-triangular entries s^, 1 < i < j < n are i.i.d. random variables such that s^j = 1 
with probability ^p, Sjj = — 1 with probability ^p, and s,j = with probability 1 — p. 

It is well known that all the eigenvalues of S n are purely imaginary numbers. Assume that 
i\i,i\2, ■ ■ ■ ,i\ n are all the eigenvalues of S n where every is a real number and i is the 
imaginary unit. Let S' n = (—i)S n . Then S' n is an Hermitian matrix with eigenvalues exactly 
Ai, A2, . . ■ , A n . Therefore, the skew energy £s(G a (n,p)) can be evaluated once the spectral 
distribution of the random Hermitian matrix S' n is known. 

Usually, it is more convenience to study the normalized matrix M n = ^pS*^ = [my]. Ap- 
parently, M n is still an Hermitian matrix in which the diagonal entries ma = and the upper- 
triangular entries m^-, 1 < i < j < n are i.i.d. copies of random variable £ which takes value 
with probability ^p, — with probability ^p, and with probability 1 — p. It can be verified 
that the random variable £ has mean 0, variance 1, and expectation 



if s is odd; 

ifs = mod 4; (3.1) 



1 



if s = 2 mod 4. 



Observe that if p = o(l), then the matrix M n is not a Wigner matrix since the moment 
is unbounded as n — > 00, and thus the limiting spectral distribution of M n cannot be directly 
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derived by the Wigner's semicircle law. However, by the moment method, we can establish 
that the empirical spectral distribution of -^=M n also converges to the standard semicircle 
distribution, which in fact generalize the Wigner's semicircle law to a larger extent. 



Theorem 3.1 For p = £*;(-), the empirical spectral distribution (ESD) of the matrix -^M n 
converges in distribution to the standard semicircle distribution which has a density p sc (x) with 
support on [—2, 2], 

Psc(x) ■= + 



x 2 . 



We first estimate the skew energy £s(G a (n,p)) by applying the theorem above but leave the 
proof of the theorem at the end of this section. Clearly, —t=Xi, -77A2, . . . , and -^=Ai, -4=A2, 

yP yP yP y n P y n P 

. . —)=\ n are the eigenvalues of M n and -4=M n , respectively. By Theorem 13.11 we can deduce 

v n P v n 



that 



£ s (G°(n,p)) 



^3/2^1/2 j^3/2jjl/2 Z— 1 

i=l 



1 n 



1 U 1 

-Y\— 

n ' . fir 



^np 



\x\dF%~ 1/2M "(x) 



2vr 

8 

3^ 



x\ p sc (x) dx (n — v 00) 

— / Iar|v4~— x 2 dx 
!tt J-2 



Hence, we can immediately conclude that 



Theorem 3.2 For p = co(-), the skew energy £s{G a (n,p)) of the random oriented graph 
G a (n,p) enjoys a.s. the following equation: 



£ s (G>,p)) = n 3 /y/ 2 (I; + o(l)) 



Now we are ready to prove Theorem 13.11 

Proof of Theorem\£j]; Let W n = -±=M n . To prove that the empirical spectral distribution 
of W n converges in distribution to the standard semicircle distribution, it suffices to show that 
the moments of the empirical spectral distribution converge almost surely to the corresponding 
moments of the semicircle distribution. 

For a positive integer k, the k-th moment of the ESD of the matrix W n is 

J x k dF^(x) = (Trace(H^) 
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and the fe-th moment of the standard semicircle distribution is 

x k p sc (x)dx. 

2 

Hence, we need to prove for every fixed integer k, 

r-2 



-E (Trace(W^)) — > j x k p sc (x)dx, as n 



oo. 



On one hand, we can determine that 

for k = 2m + 1, f\ x k p sc (x)dx = due to symmetry; 

for k = 2m, 

/2 1 r 2 i r 2 

x k p sc {x)dx = — / x k y / 4 - x 2 dx = - x 2m \A - x 2 dx 
2 2?r J-2 n Jo 

n2m+l i-l 

= / y m ~ 1/2 (l - y) l/2 dy (by setting x = 2jy) 

2 2m+i r(m + l/2)r(3/2) 1 (2m 



it T(m + 2) m + l\m 

On the other hand, we expand the trace of W k into 

iE (Trace(^)) = -^E (Trace^ 



J— ^ E(7n ili3 mi a i 3 ---m ifcil ). (3.2) 



n l+fc/2 

!<«!,. ..,ife<n 



Every term in the sum above corresponds to a closed walk of length k in the complete graph 
of order n. Recall that the matrix M n satisfies that the entries m™, 1 < i < j < n, are i.i.d. 
copies of the random variable £, which commits |E(^ S )| = if s is odd and |E(£ S )| = ^^-^ s _^ if 
s is even. Besides, my = mj7 = —rriji. For convenience, we also regard mjj as an edge and m^, 
the inverse edge of m^-, or vice versa. 



When k is odd, each walk in the Sum (13.2p contains such an edge that the total number of 
times that this edge and its inverse edge appear in this walk is odd. Apparently, by Equ. (|3.ip 
and the independence of the variables, this term is zero. Thus 

-E (Trace(W^)) = 0. 

When k is even, suppose k = 2m and let t be the number of distinct vertices in a closed 
walk. All closed walks in the Sum (|3.2|) can be classified into the following two categories: 

Category 1: There exists such an edge in the closed walk that the total number of times that 
this edge and its inverse edge appear is odd. Similarly, by Equ. (|3.ip this term is zero. 



Category 2: Each edge in the closed walk satisfies that the total number of times that this 
edge and its inverse edge appear is even. It is clear that the number of distinct vertices in this 
walk t < m + 1. We then continue to divide those walks into the following two subcategories: 
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Subcategory 2.1: t < m. It is clear that the number of such closed walks is at most n t ■ t k . 
Then these terms will contribute 



- m 

t=i \{i 1 ,...,i k y\=t 

2m-2(t-l) 



^ 1 m t ' ^ 

<— - — - m- n ■ m 

- n i+m 

_m t+1 _ ( 1 



2m-2(m-l) 



The first inequality is obtained by merging the same edges and their inverse edges together and 
then employing Equ. (|3.ip . The second inequality is due to the monotonicity. 

Subcategory 2.2: t = m + 1. In this case, each edge in the closed walk appears only once, and 
so does its inverse edge. By E(££) = — E(£ 2 ) = 1 and the independence of the variables, this 
term is 1. And the number of such closed walk is given by the following lemma. 



Lemma 3.3 \2j The number of the closed walks of length 2m which satisfy that each edge and 

-J- l 

rn+l 



its inverse edge in the closed walk both appear once is ^Vj ( 2 ^ 1 ) . 



m+l \ m ) \ np 



Prom the above discussion, it follows that 

1 / u \ fo if fc = 2m + l; 

-E(Trace(^ n fc )) = { „ / \ 

which implies that if p = w(-), then 

-E ^Trace(P^)^ — > J x h p sc (x)dx, as n — > oo. 

The proof is thus completed. 



4 The skew energy of G a nd 

In this section, we consider the skew energy of random regular oriented graphs. We first 
recall [3] the definition of a random regular graph G n ^, where d = d(n) denotes the degree. 
G ni d is a random graph chosen uniformly from the set of all simple (f-regular graphs on n 
vertices. A random regular oriented graph, denoted by G a n d , is obtained by orienting each edge 
of the random regular graph G n> d, randomly and independently, with probability 1/2. Let A n 
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be the adjacency matrix of G n ^d and R n be the skew-adjacency matrix of G a n d . The estimates 
of the skew energy of G° d are different in the cases of d fixed and d — > oo. Therefore, we shall 
discuss these two cases separately. 

4.1 The case that d > 2 is a fixed integer 

In this subsection, we estimate the skew energy of G^ d , where d > 2 is a fixed integer. We 
first recall the fact about the limiting spectral distribution of the random regular graph G n< d 
with d fixed (which means the limiting spectral distribution of the adjacent matrix A n ), which 
was derived by McKay [13J . 

Lemma 4.1 ]13^ Let G n) d be a random regular graph with the adjacency matrix A n . If the 
degree d is a fixed integer and d > 2, then the empirical spectral distribution F„ approaches the 
distribution F(x) whose density function is 



Remark 4.1. McKay [13] used the moment method to prove the lemma above, i.e., he proved 
that for each k, the fc-th moment of the ESD of the matrix A n converges to the fc-th moment of 
the distribution F(x), 



Note that Trace(A^) is the number of closed walks of length k in A n . When d is fixed, the graph 
G n ^ is almost surely a locally (i-regular tree. 

Now we consider the random regular oriented graph G a n d . Set T n = {—i)R n = [Uj]. For a 
fixed A;, the limit of the fc-th moment of the ESD of T n is = lim (Trace(T^)) . We note 
that when d is fixed and n — > oo, the underlying graph is almost surely a locally d-regular tree. 
If one oriented edge appears in a closed walk of length k, then its inverse oriented edge appears 
with the same number of times. We can get that mk is equal to the number of closed walks of 
length A: in a d-regular tree starting at the root. Combining with Lemma 14.11 we conclude the 
following theorem. 

Theorem 4.2 Let G^ d be a random regular oriented graph with the adjacency matrix R n , and 
let T n = {—i)Rn = [tij]. If the degree d is a fixed integer and d > 2, then the empirical spectral 
distribution F^ n approaches the distribution F{x) which has the density function pd with support 





on 



[-2 v / d^T,2 



Vd^T] 



Pd = 



d^/A{d - 1) - 
27r(d 2 - x 2 ) 
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We now turn to the estimate of the skew energy £s(G^ d ). Note that £s(G^ d ) also equals the 
sum of the absolute values of all the eigenvalues of R n . Suppose Aj., A2, • • • , X n are the eigenvalues 
of Rn. By Theorem 14.21 we can deduce that 



n n 

i=l 



1 n 

n ^— ' 

= 1 

\x\dF^ n {x) 

\x\ pd{x) dx (n — >• 00) 

Z^ 2 ^ 1 <y4(rf- 1) -x 2 

J 2ir(d 2 - x 2 ) 

2dVd^T d(d - 2) 2\/d 3 T 

= • arctan — . 

7T 7r a — 2 

To summarize, we can obtain the following theorem. 

Theorem 4.3 For any fixed integer d > 2, the skew energy £s{G a nd ) of the random regular 
oriented graph G a d enjoys a.s. the following equations: 



„ . f2dVd^T did -2) 2^/d^l 
£ S {G a n4 ) = n —± ^ >- ■ arctan y + o 1 

\ TT TT CI — 2 

In particular, when d = 2, £s(G^ d ) = n (A/it + o(l)). 



4.2 The case that d — > 00 



The estimate of the skew energy of G^ d with d — > 00 depends on the following key lemmas. 

Lemma 4.4 \17f If np —> 00, then the random graph G(n,p) is np-regular with probability at 
least exp (— 0(n(np) 1//2 )) . 

Next, we consider random oriented graphs. By the definitions of a random oriented graph 
and a random regular oriented graph, we can generalize the lemma above into a result for random 
oriented graphs as follows. 

Lemma 4.5 If np — > 00, then the random oriented graph G a (n,p) is np-regular with probability 
at least exp (— 0(n(np) 1 / 2 )) . 

Lemma 4.6 Ji7| / Let M be an n x n Hermitian random matrix whose off-diagonal entries 
£ij are i.i.d. random variables with mean zero, unit variance and |£y| < K for some common 
constant K. Fix 5 > and assume that the fourth moment M4 := sup^-E(|£y| ) = o(n). Then 
for any interval I C [—2,2] whose length is at least Q(5~ 2 / 3 (Af^n) 1 / 3 ) , there is a constant c 
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such that the number Nj of the eigenvalues of which belong to I satisfies the following 

concentration inequality 



N z -n I psc(t)dt 



> 5n J p sc (t)dt\ < 4exp ^— c 



K 2 



Consider the random oriented graph G a (n,p) with np — > oo as n — > oo and the skew- 
adjacency matrix S n . Recall that M n = ^S n . For an interval / let N'j be the number of 
eigenvalues of M n in /. Apparently, M n satisfies the condition of Lemma 14.61 (M := M n , K := 
1/ yjp). Thus, we can immediately obtain the following lemma. 

Lemma 4.7 For any interval I C [—2,2] with length at least ( ^4°^i/2 ) 1//5 ; we have 



x 



N'j-n J p sc (x)d. 
with probability at most exp(— cn(np) 1 / 2 log(np)) 



> 5n 



J p sc (x)da 



By Lemmas 14.51 and 14.71 the probability that N'j fails to be close to the expected value in 
the model G a {n,p) is much smaller than the probability that G a (n,p) is np-regular. Thus, the 
probability that N'j fails to be close to the expected value in the model G a n d where d = np is the 
ratio of the two former probabilities, which is 0(exp(— criy/nplognp)) for some small positive 



constant c. Recall that R n is the skew- adjacency matrix of Set L r 

be the number of eigenvalues of L n in /. Thus we can conclude that 



/ d/n 



-Rr, and let N'J 



Theorem 4.8 ( Concentration for ESD of G^ d ) Let 5 > and consider the random regular 
oriented graph G a nd . If d tends to oo as n —> oo, then for any interval I C [—2,2] with length at 
least S- 4 ^d- 1 / 10 loR 1/5 d, we have 



N'j - n / p sc (x)dx 



< 5n J p sc (x)dx 
with probability at least 1 — (9(exp(— cny/d\og(d))) . 

Theorem 14.81 immediately implies a result as follows. 

Theorem 4.9 If d -too, then the ESD ofrT 1 ! 2 ! n converges to the standard semicircle distri- 
bution. 



Now we are ready to estimate the skew energy of the random regular oriented graph G° d . 
Suppose that 1X2, • • • , i\ n are the eigenvalues of R n . Then, (f -1 / 2 Ai, (i~ 1 / 2 A2, . . . , d~ l l 2 \ n 
are the all eigenvalues of n~ l / 2 L n . By Theorem 14.91 we can deduce that 
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a.s. 




x\ p sc dx (n — > oo) 



x|v4 — x 2 dx 
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Therefore, the skew energy £s(G^ d ) can be formulated as 



e s {G° td ) = nd 




We can thus immediately obtain the following theorem. 

Theorem 4.10 For d = d(n) — > oo, the skew energy £s(G° d ) of the random oriented graph 
G a nd enjoys a.s. the following equation: 
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